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Abstract. The increasing interest in nonequilibrium effects in condensed matter theory 
motivates the adaption of diverse equilibrium techniques to Keldysh formalism. For methods 
based on multi-particle Green or vertex functions this involves a detailed knowledge of 
the real-time properties of those functions. In this review, we derive general properties of 
fermionic and bosonic multi-particle Green and vertex functions for a stationary state described 
within Keldysh formalism. Special emphasis is put on the analytic properties associated 
with causality and on a detailed discussion of the Kubo-Martin-Schwinger conditions which 
characterize thermal equilibrium. Finally we describe how diagrammatic approximations and 
approximations within the functional renormalization group approach respect these properties. 
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1. Introduction 

Multi-particle Green and vertex functions are an important tool for the theoretical analysis 
of interacting quantum systems in condensed matter physics, n-particle Green functions 
contain the relevant information to compute n-particle properties of the system 1 1 1. Moreover, 
determining multi-particle functions may be an auxiliary step on the way to find the single- 
particle Green function (or an approximation to it) and single-particle properties. For example, 
the functional renormalization group (fRG) approach [2] leads to a system of flow equations 
which couples all «-particle functions with n > 1 in an infinite hierarchy. Depending on 
the particular fRG scheme, these n-particle functions may be Green or irreducible vertex 
functions. Another example is given by diagrammatic resummation techniques based on 
the Bethe-Salpeter equations, such as fluctuation exchange |3| or parquet-equations [4, 5|. 
For these approaches the preliminary task is to determine the irreducible two-particle vertex 
function, from which in turn the single-particle self-energy can be deduced. 

The methods mentioned above are typically used within zero-temperature or equilibrium 
formalism. Calculations in equilibrium formalism are based on imaginary frequencies and 
require an analytic continuation to real frequencies in order to determine many physical 
observables of interest. At finite temperature, such analytic continuation is numerically quite 
nontrivial even for a single frequency (see e.g. IS] and references therein). Therefore a 
formulation of the problem within the real-time (real-frequency) Keldysh formalism |7 8 9| 
may be advantageous. Beyond that, Keldysh formalism allows for the description of 
nonequilibrium situations such as stationary state transport at finite bias voltage. 

A central ingredient of Keldysh formalism is a time contour consisting of two branches, 
one of them describing evolution forward and the other backward in time. As a consequence of 
the existence of these two branches, the n-particle Green function is a tensor in 2-dimensional 
space of rank 2n, having 2^" components, each being a function of 2n state indices and 2n 
frequencies [or (2n — 1 ) frequencies in the case of time translational invariance] . The resulting 
additional complexity of the formalism is however reduced due to the fact that the components 
are not independent from each other but obey certain relations connecting them. In this 
review, we give a comprehensive and self-contained presentation of fundamental relations 
for the multi-particle Green and vertex functions within Keldysh formalism: the behaviour 
under exchange of particles and under complex conjugation, properties connected to causality 
and the Kubo-Martin-Schwinger (KMS) conditions characterizing thermal equilibrium. The 
special single-particle versions of these relations are commonly known, e.g. the fact that 
retarded and advanced single-particle Green function are adjoint to each other, and the 
fluctuation dissipation theorem in thermal equilibrium. However, also the general multi- 
particle versions are required to adapt the above-mentioned methods to Keldysh formalism, 
for instance for real time implementations of the fRG, which have been initiated in ifTOlfTTl . 

Some of the properties in question have been analysed in lfT2l[T3]| for multi-point Green 
functions based on real bosonic fields. The adaption of their results to condensed matter 
problems is far from being obvious due to two major complications. Firstly, condensed 
matter problems usually involve complex bosonic or Grassmann fermionic operators instead 
of real bosonic ones. This makes necessary a separate book-keeping of incoming and of 
outgoing indices (corresponding to annihilation and creation operators) instead of just one 
type of indices as for real bosons. Additionally, due to the commutation relations of complex 
bosonic or Grassmann fermionic operators, the initial order of a time-ordered operator product 
is relevant for the result, as opposed to real bosons. Furthermore, the time reversal properties 
of the field operators which are important for the analysis of the KMS conditions differ 
from those of real bosonic field operators. The second source of complication is due to the 
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fact that condensed matter problems usually involve external potentials which lift diverse 
symmetries of the system. Reference ifTJI implicitly assumes translational invariance in space 
and time and invariance under the product of time reversal and parity. These prerequisites, 
which considerably simplify the investigation of the KMS conditions, are not given in typical 
condensed matter problems. Despite these fundamental differences some of the concepts 
developed in references lfT2l[T3]| can be fruitfully adapted to our analysis. 

In sectionlH we start by a short description of the physical situation under consideration. 
The set-up of an interacting quantum system coupled to a number of noninteracting reservoirs 
represents a typical stationary state configuration described in Keldysh formalism. 

We continue by specifying precisely the definitions of the Green and vertex functions 
in question in section |3] Here it is necessary to take care of the details (like prefactor 
conventions) since the exact form of the relations to be discussed depends on them. 

In sections 14.11 and 14.21 we derive useful relations connected to the permutation of 
particles and to complex conjugation. Section 14.31 focuses on a theorem of causaUty lfT2l 
which we cast in concise form. We show that this theorem has direct consequences for the 
analytic properties of the Green and vertex functions in Fourier space and identify components 
which are analytic in one half plane of each of their (2n — 1 ) independent frequency 
arguments. Remarkably, our derivation does not use the Lehmann representation. This fact is 
highly advantageous due to the unwieldy form the latter acquires in the multiparticle case. 

Section |5] is devoted to the special situation of thermal equilibrium. The fact that the 
density matrix then corresponds to a time evolution in imaginary time is at the root of the 
KMS conditions |fT4l[T5|. In section ISTl we describe how to concisely formulate the multi- 
particle KMS conditions exploiting the so-called tilde Green function which is defined on the 
analogy of |13| by means of a non-standard ordering of operators on the Keldysh contour. 
However, in order to extract useful information from the KMS conditions, the tilde Green 
function has to be eliminated in favour of the standard Green function. In sections |53] and l534l 
we find that time reversal can be used for this purpose, an observation similar to the one 
of lfT2l for multi -point Green functions based on real bosons. In section |53] we characterize 
an important class of systems (including e.g. the Anderson impurity model 1 16 |) with specific 
time reversal properties. For this class, we formulate a multi-particle FDT in section lS^l 

Finally we study under which conditions approximations to the Green or vertex functions 
conserve the before-mentioned properties (section |6]l. These conditions turn out to be quite 
handy for the large class of diagrammatic approximations. For approximations within the 
fRG framework, the choice of the flow parameter is found to be crucial for the conservation 
of causal features and of the KMS conditions. 



2, The system under consideration 

The physical setup underlying our considerations is a finite interacting quantum system of 
a single kind of either bosonic or fermionic particles which is coupled to M noninteracting 
reservoirs, M > 1. Let {q] — {s} U {ri} U . . . U {^m} be a basis of single-particle states which 
are either localized in the system (states s) or in one of the reservoirs (states r^, k — 1, . . . ,M). 
We use the index q to refer to any of those states. The total Hamiltonian under consideration 
is given by 

H — Hsys + 52 ^res + 52 ^coup, (1) 

k=l k=\ 
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with 



4 ^"i-<2l*l^'2 
12,4 



Hml = h^^\^^a^^ar^, (3) 

^/iolp = L''.vk,«K+H.c. (4) 

in standard notation of second quantitation. Here 

hq'\c,^ {q'\h\q) = h*\y, (5) 
are the matrix elements of the single-particle Hamiltonian and 

Vs[s'2\hs2 = (sWiMsisi) + l^{s[s'2\v\s2Si) (6) 
denote the (anti-)symmetrized matrix elements of the two-particle interaction, with 

{+1 for bosons, 
— 1 for fermions. 

At some initial time fo the whole configuration is assumed to be described by the density 
matrix 

P (fo) = Psys <E) Prei «)...«) p4f ' , (8) 

where 

PkI ^ WTO ^) ^=1,...,M, (9) 

is the grand-canonical distribution of reservoir k with temperature 7). = 1 / and chemical 
potential jj.^. Keldysh formalism allows us to compute properties described by the density 
matrix p(f), t > to, which evolves unitarily from p(fo) under the full Hamiltonian H given 
iniBlIZlSia. 

In the limit of an infinite reservoir size L and an infinite transient time t — to such that 
L/u^ —to — > °° (with u being the particle velocity in the reservoir), expectation values under 
p{t) become stationary and do not depend on the initial density matrix Psys of the finite 
interacting region. We will restrict our considerations to this stationary situation. 

A particular stationary situation is a global equilibrium which occurs if all and jj.). 
are equal. We specialism to the equilibrium case in section |5] The results found there are 
also applicable to interacting bulk systems (that is large systems with the Hamiltonian (|2]l) in 
thermal equilibrium. 



3. Definition of tlie Green and vertex functions 



The time- or frequency-dependent n-particle Green function is defined by 

Gl^^^Mt')^Hy\Tj,/,'\h)---al/:Ht„^^^^^^ (10) 

Gjj{,(a)|a)') = |dfi...df:ei('»-'-"'-'')Gjj{,(f|f'), dD 

where q—{q\,.. .,qn), f = (fi , . . .,?„), (0 — {(Oi , . . . , »„) and j ~ {ji ,...,]„) are multi-indices 
denoting state, time, frequency and branch of the time contour, respectively. The annihilation 



Properties of multi-particle Green and vertex functions within Keldysh formalism 



5 



A\ A2 A3 




+ • • • — -' t 

Aa As Afi 



Figure 1. The time contour consisting of the — and the +-branch with some example 
operators Ai , . . . ,A(,. The time ordering operator 7^ will establish the order A4A5A5A3A2A1 . 
The operator used in the definition )411 of the tilde Green function will establish the order 
A6A5A4A1A2A3. 

and creation operators and aj^ of a particle in the state qj^ obey standard commutation 
relations 

[aq^^a'li^_j^ — aq^a^^ — ll,a^ij^aq^ ~ 8q^^q2^ (12) 
[aq,,aq,]_t^ = [a\^,a\^]_l^^Q, (13) 

and aqf\tii) denotes an annihilation operator at time tj^ in the Heisenberg picture with 
reference time to. The operator is situated on the branch j/^ = =F of the time contour |7|. 
The contour ordering operator in (fTOl i rearranges the sequence of operators in the following 
way: all operators with +-index are placed left of all operators with — index. The block of 
operators with +-index is anti-time ordered internally (which means time arguments of the 
+-operators increase from left to right), whereas the block of operators with — index is time 
ordered (which means time arguments of the — operators decrease from left to right). By 
definition of 7^ the expression ordered in this way has to be multiplied by a factor ^, if it 
represents an odd permutation of the initial order given in ( fTol i. The branches of the contour 
are depicted in figure [T] where also an example for the order established by Tq is given. The 
expectation value in ( fTOl l is defined by 

(A(f))=Trp(fo)A(f), (14) 

with p (to) being the density operator at time to^ In (fTTT i we used the shorthand notation 
(O t = (Oiti + . . . + 0}„t„. 

Due to the time translational invariance of the stationary state, the «-particle Green 
function can be expressed as a function of only (2n — 1 ) independent time or frequency 
arguments. For example f 1 and (0\ can be eliminated by 

G{ti,...,t„\t[,...,t'„)=G(0,t2-tu...,t„-ti\t[ -fi,...,f;,-fi), (15) 

G((Oi, . . . ,(o„|(oJ, . . . , «,') = 2n5{coi + ... + «„- «J - ... - o),') 

xG(ri ^0,co2,...,co„\co[,...,co'„). (16) 

However, many of the properties discussed in the following are formulated most easily when 
all time and frequency arguments are treated on equal footing. Therefore we will use the 
reduced number of frequency arguments only occasionally, e.g. for the study of analytic 
properties in section |43] 

Apart from Green functions we will also discuss the properties of one-particle irreducible 
vertex functions. These vertex functions can be derived from a generating functional called 
effective action which is the Legendre transform of the generating functional for the connected 
Green functions. Details on this path-integral approach to the irreducible vertex functions can 
be found for Matsubara formalism e.g. in 1 1] and for Keldysh formalism in 1 10|. Since we do 
not need any recourse to the generating functionals in the following, we choose an equivalent 
way to introduce the vertex functions which is based on a diagrammatic expansion. 
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To set up a standard Hugenholtz diagrammatic technique Q] for a problem with 
instantaneous two-particle interaction v, we define the (anti-)symmetrized interaction vertex 

or, in Fourier space, 

with v^'j^^|£,,^2 given by (|6]l. Since a single Hugenholtz vertex incorporates both the direct and 
exchange Feynman vertices, this diagrammatic language is especially concise in formulation. 
According to Wick's theorem, the n-particle Green function can be expanded into a sum of 
diagrams which consist of lines representing the noninteracting one-particle Green function 
G^''' and of interaction vertices ifTTl . 

For the discussion of the vertex function we will restrict ourselves to the frequency 
representation. The (anti-)symmetrized irreducible n-particle vertex function 

r%i^\^-4"^:;'^^^ (19) 

can be defined diagrammatically as the sum of all one-particle irreducible diagrams with 
n amputated incoming lines (having state, frequency and contour indices q,CO,j) and n 
amputated outgoing lines (having indices q' ,(£)' ,f). In order to evaluate a given diagram, 
one determines first the symmetry factor S which is defined as the number of permutations 
of vertices mapping the diagram onto itself. Then one chooses arbitrarily at each vertex the 
order in which the attached incoming (outgoing) lines are assigned to the incoming (outgoing) 
indices of the vertex. This order should be kept fixed during the further evaluation. Let 
count the pairs of equivalent lines in the diagram, where two lines are called equivalent if 
they connect the same vertices and run in the same direction. One can define «ioop to be the 
number of internal loops made up from internal propagators. (Inside a given vertex Vi/2'|i2 
we assume that one line continues from 1 to 1' and the other continues from 2 to 2'.) Finally, 
the incoming line k being connected to the outgoing line P(fe') via internal lines and vertices 
determines a permutation P of { 1 «'}■ The value of the diagram is then given by 



2'''=')5 



;n-l 



np^F-jnc'"'. (20) 

where one has to integrate over all internal frequencies and sum over all internal states and 
contour indices. Due to the (anti-)symmetrization of the vertex in (fTSl ). the overall sign of the 
diagram is actually independent of the chosen assignment of lines to indices at each vertex. 
An example for the application of the diagram rule is given in figure |2] The Hartree-Fock 
contribution to the one-particle vertex depicted in figure |2](d) for instance is evaluated as 

C(2^)'7i I L j Am^O^-v^'^;^'^^^^^^^^^ (21) 

Given time translational invariance one often works with vertex functions depending on 
(2n— 1) frequencies only, e.g. y{t[ = Q , (1)2 ■ ■ ■ (£>'„\(i)i ...a)„). The full vertex function ( fT9] l is 
then given by 

y{(0[ ,...,(0'„\(0i,...,(0„)= 2Kd{03[ + ... + (0'„~(0]_- ...-(On) 

y.y{t[ =0,a)^,...,a),^|coi,...,co„). (22) 

The diagrams for y{t[ ^ Q , 0)2 ■ ■ ■ (£>'„\(£i\ . . . (£>„) typically consist of single-particle Green 
functions depending on a single frequency [occurring in the second line of (fT6] l1 and 
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Figure 2. Examples for the application of the diagram rules. The lines are free one-particle 
Green functions, the dots are interaction vertices. The indices represent contour index, state 
and frequency, e.g. 1 = {ji,q\,(Oi). Diagram (a) contributes to the two-pai1icle vertex 
function. It has one pair of equivalent lines, n^q = 1 , and no symmetries, S = 1 . (b) and 
(c) show two possible, equivalent ways to assign the lines to the indices of the vertices. In (b) 
the vertex at the bottom e.g. is Vi'3'|i3, and nioop = 2, f = 1. In (c) the vertex at the bottom is 
^l'3'|3i = ?*'i'3'|l3> 'Jloop = 0, f = f . Compared to (b), three of the four vertices acquire 
a factor ij; taking into account the prefactor ^"loop^'' the possibilities (b) and (c) hence lead 
to the same value of the diagram. The Hartree-Fock diagram (d) (iieq = 0, 5 = 1) contributes 
to the one-particle vertex, that is the self-energy. The figure shows the evaluation with vertex 
vy2'\n andnioop = 1- 



of frequency independent vertices [the part {~ ji)^}i=p=j[=i'^'^q\q'^\mq2 of (flSTll. When 
evaluating such diagrams, the factor [2nY" in ( |20] | has to be replaced by (27r)""' and i/ {2n)^ 
by i/ {2n). For the Hartree-Fock diagram in figure|2](d) e.g. this yields 

C— E E /d«2/;1'/'''^GWfli^(co2). (23) 

Equation (l20b defines the prefactor of the vertex functions in such a way, that they can be 
used themselves as vertices in diagrams with the identical prefactor rules being applicable as 
for the bare n-particle interaction vertices. The self-energy is equal to the one-particle vertex, 
S = 7«=i- 

In the one-particle case it is common to call the components of the Green and vertex 
function chronological, lesser, greater and antichronological, 

G'= = G"I", G<=G"I+, G>=G+'", G'^ = G+I+, (24) 
E'= = 7-l-, E< =7"l+, E> ==7+'", E'^ = 7+l+, (25) 

respectively. 

Apart from the contour basis with indices ji^ = =F we will use the Keldysh basis ifTSll with 
indices ak = 1,2 whose components have a direct physical interpretation in terms of response 
and correlations ||T9l . The transformation to the Keldysh basis is defined by a change of 
operators, 

^(D^fl a ^(2) ^ « +« ^ (26) 

V2 V2 
The Green and vertex functions change according to the tensor transformation laws, 

G«l«' = (D-1)"I^G^'I^'d-'"I"', (27) 
7«'l« = (D-')«'I^'7^"'^''D^''", (28) 
where an impUcit summation over all internal indices is assumed, and 

D^I« = flD«l«*, (29) 
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(D-')'l- = (D-'fl^ = -lj, (D-')'l+ = --L. (31) 

In the particular case of one-particle functions, we have 

G^l^=0, g'|2 =G^\ G^^^=G^'\ G^\^^G^, (32) 

/|i=E'^, /|2=E^^\ 7211 =eA\ y2|2^Q^ (33) 

where the retarded, advanced and Keldysh component of the Green function and the self- 
energy are given by 

G^^'(f|O--i0(f-O(k(O,«;(O]_f), (34) 

Gtl, (f |f ') = i0(f ' - ( k (0 , 4 (f ')] _^ ) , (35) 



G 



and 



(jRet ^ ^(O)Ret ^ ^(O)Ret ^ j^Ret „ ^Ret^ (37) 

= g(°'^^ + g'")'^" o E^" o G^^ , (38) 
=g'^'='oE'^oG^\ (39) 

Here we use "o" as a shorthand for summation over internal states and integration over internal 
times. 



(AoB)^l^,(f|f') / dsA^\Ms)Brw{s\t'). (40) 

r •' 

When discussing the KMS condition in section 15.11 we will need yet another type of 
Green function, the so-called tilde Green function which we define on the analogy of lil3J . It 
is given by 

&l(^m = i-^T{fJ^!\h)...a^^:\tn)af^ (41) 

G%((0\(0') - /df,...df:e'(-'-»'-'')G;l|^(f|0, (42) 

which differs from ( fTOl l and (fTTI ) only by the other ordering operator 7^. Like Tc, 7^ places 
all operators with +-index left of all operators with — index. However, under 7^ the block 
of operators with +-index is time ordered internally, whereas the block of operators with — 
index is anti-time ordered. Again, if the final order is an odd permutation of the initial one 
it has to be multiplied by a factor C,. An example for the order established by Tq is given in 
figure [U 

For the one-particle case it is straightforward to check the identity 

G%{t\t')^Glll,W). « = 1, (43) 

where a bar over a contour index indicates a swap of the branch, j = —j for 7 = T- However, in 
the general multi -particle case there is no direct possibility to express the tilde Green function 
through an ordinary Green function since the orders established by and 7^ are inherently 
different; we will need to make use of time reversal to connect both quantities, see section |53] 
A direct consequence of having a special relation for n = 1 is the particular structure of the 
single-particle FDT, see section]; 
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4. General properties of the Green and vertex functions 

This section is devoted to properties of the Green and vertex functions which are valid in a 
general steady state, whereas section|5]focuses on the particular case of thermal equilibrium. 
We will derive the properties of the Green functions starting out from ( fTOb and ( fTTT i. Once we 
know a certain relation for the Green function we can study its implications for the amputated 
irreducible diagrams, the values of which are given by ( |20] i. This will allow us to infer a 
corresponding relation for the irreducible vertex functions. 

4.1. Permutation of particles 

Let P be a permutation of (!,...,«) and for any multi-index b = {bi,...,b„) define Pb : = 
{bp(j),...,bp^„-j). The prefactors ^ due to time ordering in dTol ) and due to permutations 
in (I20I 1 entail 

g';%,{P(o\(o') = G%{,{co\Pco') = rGjl{,(«|to'), (45) 

7S'|^(Pa)» = T^:|;>'|Pa)) = r^^l^'lco). (46) 

The identical equations hold after transforming to the Keldysh basis. The contour indices j 
are merely to be replaced by Keldysh indices a. 

4.2. Complex conjugation 

Conjugating (fTOl l interchanges creation and annihilation operators and reverses the order of 
operators. The reversed order is identical to the order obtained by Tc when each contour index 
iii is swapped to its opposite value —ji,. Therefore, we find 

G%ityr = {-\fG^;)l{t'\t) (47) 

with j — {ji,-.. in) and ji^ ^ —j^. Note that the parity of the permutation needed to establish 
the order according to Tc is identical for G'''/, (r|f') and G'' j/ {t'\t). Hence no additional sign 
prefactors are required in i47i . Fourier transformation leads to 

G^^,{co\cor^{~irG^}l{co'\co). (48) 

In order to derive the corresponding identity for the vertex function we discuss how 
to conjugate the individual diagrams which the vertex function is composed of and whose 
values are given by ( |20] |. The noninteracting single-particle Green functions entering the 
expression ( l20b can be conjugated according to ( l48T l. The conjugation rule for the bare two- 
particle interaction vertices in expression (l20l i is inferred from their definition (fTSl ) to be 

= il^'i<o'\<o) = vi|;:(«l«') = -y%{o,\o/), (49) 

where we make use of {q\q'2\v\q2qi) — (?2?i I v|^i^2) in the first step. Applying ( |48] | and ( l49b 
to the conjugate of (l20b we find that complex conjugation maps one-to-one the diagrams 

contributing to 7j,|'|^(to'|(o) onto those of T^j^, {co\(o'). The precise relationship can be deduced 
from the fact that an amputated n-particle diagram comprising m two-particle vertices contains 
(2m — n) free propagators, and is found to be 

r!l^{(o'\(or^-r!^{(o\(o'). (50) 
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When transforming to the Keldysh basis via ( |27] l and (|28] l we use that D is real and fulfils 
the relation 

= (51) 

which follows from ( |29] l - (ISTT l. We find 

G^";l«' = (-l)"+^*(«*+<)G;;|l"(to'|«), (52) 

7;;i"(«»* = (-i)^+£*(«*+«*)7;;i?Vi«'). (53) 

For the one-particle functions, equations (l52b and (|53] | yield the well-known relations 



G^'^,{(o\(o') = gX{(o'\(o)\ (54) 



^%{co\co') = G% 
G^^,(co|co') = -G^I^(«»*, (55) 

E^^;,(«|«')=E^;(«'|co)*, (56) 
I^\^,{(o\(o') = -i:%{(o'W. (57) 

4.5. Causality and analyticity 

Consider the Green function G-'l-''(f |f') with fixed sets of times f = (fi , . . . ,f„), f' = (fj, . . .f^), 
where fj happens to be the greatest time argument, fi > f2, . . . ,f^. In this case the order of 
operators established by Tq is independent of the contour index ji: given ji — — , a^' (fi) is 
sorted to the leftmost position of all operators with — index; given ji = +, a-" (fi) is sorted 
to the rightmost position of all operators with +-index which is exactly the same place in the 
total expression as for ji = — . As a consequence 

G-J2...7„l/(f|,')^G+J2...;„l/(,|f')^ if ti>t2,...,t'„- (58) 

This property is at the heart of a theorem of causality which emerges after transforming ( fSST l 
to the Keldysh basis. Let us first transform only the contour index ji to the Keldysh basis: 

Qa,.j2-j„\j\t\t') = (D-i)"il-G-'^2-^"l^"(f|f') + (£>"')"'l+G+'^2-^"l^"(f|f') 

= [(D-')«'l- + (D-')«il+]G-'^2-^"l^-'(r|f'), if fi >f2,. (59) 

For ai = 1 we can use 

(D-')i|---L = -(Z)-i)i|+ (60) 

to see that ( |59] ) vanishes. After transforming the remaining contour indices to the Keldysh 
basis, we find 

Gi'"2-""l«'(f|r') = 0, if ti>t2...t'„. (61) 

The same line of argument holds if any other time from ti,...,t'„ is the maximum time. [If the 
maximum time is one of fj, . . . we use D^l^ = — D+l' instead of (l60ll.l This leads us to the 
following theorem of causality lfT2l : the Green function G"l"'(f|f') vanishes if the Keldysh 
index associated with the strictly greatest time argument equals 1. (The theorem does not 
apply when multiple time arguments have simultaneously the maximum value. However, this 
particular configuration has a vanishing measure in the 2n-dimensional time space.) 
A special case is given when all Keldysh indices equal 1, 

G' ■■'I' eeO. (62) 
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If only one of all Keldysh indices equals 2 the corresponding time argument has to be 
the greatest one for a nonvanishing Green function. This has a direct impact on the analytic 
properties of the Green function in frequency space as we will show now. 

Consider the Green function (j2i...i|i . i sn{[\^ all Keldysh indices being 1 except for 
a\ = 2. Due to the theorem of causality (521...1I1. .1 (^^^ . .f„|fj ...t'^) is nonvanishing only for 
ti > f2, . . . , f/,. For the Fourier transform ( fT6] l this means that the integration range of f2, . . . 
is bounded from above by f 1 , 

^2i...i|i...i (fj ^ 0, a>2, . . . , (o„\(o') = r df2 . . .df;e'('^'2+...+'»,A,-ffl' 



xG21 ■'l'-i(fi =0,f2,...,f„|f')- (63) 

G^'' (f|f') decays if any of |fi. — fi |, — fi | goes to infinity [which makes the integral 
in (l63T l converge]. Due to the upper boundary for the integration range of f2,...,f„, only 
the decay for f2, . . . ,f„ —°° will create poles (branch-cuts) in the upper half plane (uhp) 
of a)2, - 0„. However, the Green function will be analytic in the lower half plane (Ihp) 
of these frequencies since nonanalytic features in the Ihp corresponding to decays in the 
limit t2,...,t„ — > +00 do not appear due to the restriction f2, . . . ,f„ < fi. In the same way, 
Q2i...i\i...i^f^ = 0, a>2 . . . a)„|a)( . . . to,') can be proven to be analytic in the uhp of ojj . . . a)'„. The 
differing half planes of analyticity for the co^ and (o[ stem from the different sign conventions 
for (0 and o' in the exponential e't™ '^'" ' '') of the Fourier transform. 

The same line of argument holds if any other Keldysh index from ai , . . . a,' equals 2, the 
rest being 1. Hence we can formulate a theorem of analyticity: The Green functions 

G2i-i|i-i(fi,c(}2,...,a)„i«;,... ,«;;), 

G'^' ■'(c(Ji,f2,C(j3,. . . ,co„|co(,... ,aj^), 

(64) 

G'-'\'-'^\(0u...,(0n\(0[,...Mn-2A-U<). 
Gl-l|l-12(t0i,...,t0„|t0(,..., 

are analytic in the Ihp of the frequencies (£)\ , . . . , (£)„ and analytic in the uhp of the frequencies 
(»{,..., O,' . These 2n special Green functions describe higher order response functions. Being 
formulated for bosons in time space, they are correspondingly given by the expectation value 
of a series of nested commutators [19 , 12 |. 

The irreducible vertex functions obey a theorem of causality as well: the vertex function 
7" l"(f'|f) vanishes if the Keldysh index associated with the strictly greatest time argument 
equals 2. 

The proof of this theorem is based on the structure of the diagrams contributing to 
7" l"(f'|f). Consider such a one-particle irreducible diagram t/" l"(f'|f) and suppose e.g. 
ai = 2. We show that a nonvanishing contribution to the diagram necessitates that at least 
one time argument from t2,---t'„ is greater or equal f 1 . 

The diagram is composed of bare two-particle vertices which are connected by 
internal propagators. The vertices and propagators contribute to the value of the diagram 
multiplicatively, as described in ( l20l ). Some of the vertices have one or two amputated external 
lines which correspond to the indices of the diagram. When transformed to the Keldysh basis, 
the bare two-particle vertex given in (fTTb acquires the form 

v-«l«3l«.«.(,;,,^|,,,,).^"- =^2 = r;=.D, if«, + a2 + a; + a^isodd, ^^^^ 
" else. 
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According to this equation all four time arguments of a given vertex are equal. Therefore the 
vertex carrying f i as external time argument has its other three time arguments equal to ti as 
well. If there is an external one among those three, then f i is not the strictly greatest external 
time argument. Let us now suppose that the other three time arguments are internal ones. As 
a\ = 2, at least one of the vertex' three internal Keldysh indices equals 1, see equation ( l65T l. 
There is an internal one-particle Green function attached to this Keldysh index 1 . At its other 
end this Green function is attached to a Keldysh index 2 (since G^'' =0) and a time Ti > fi; 
this follows from the theorem of causality for the Green functions. Actually all four time 
arguments of the vertex attached to the back-end of the Green function are equal to Ti > fi . 
If there is an external one among them, then ti is again not the strictly greatest external time 
argument. If not, this vertex also has as least one Keldysh index equal to 1 . Again there is 
attached a one-particle Green function leading to a vertex with time T2 > Ti >ti. We can 
proceed in this way until we arrive finally at a vertex with time T,„ > . . . > Ti > t\ which 
carries an external time argument. This external time argument T„, > f 1 concludes our proof. 
(There is the possibility that we never end at a vertex with an external line because we can get 
trapped in a loop of internal vertices with identical time arguments; this case however is of 
measure zero and vanishes when integrated over internal time arguments.) 

We mention two consequences of the theorem of causality for the vertex function which 
can be derived in the same way as their counterparts for the Green functions. First we conclude 
that 

y2...2|2...2 ^ Q (-gg-j 

Second, we get a theorem of analyticity, namely that the 2n vertex functions 

f^-^\^-\t[M2,---.<\^U---.^n). 
^12...2|2...2^^ , f2, • • • ' ' • • • ' ®«)' 

(67) 

7^- -^'^ -^'^(wj , . . . , oj/Jcoi, . . . , a)„_2,rn-i, ft),,), 

Y^-^\^-^\c0[,...,(0'„\C0u...,C0„^l,tn) 

are analytic in the Ihp of the frequencies (o[,...,(X)'„ and analytic in the uhp of the frequencies 
toi,...,©,,. 

For the special case of one-particle functions, the causality properties reduce to the well- 
known statements 

G'l'=0, I?\^ = 0, (68) 

and G^''\co') = G^^^t = 0\co') is analytic in the uhp of co', G^^ico) = G^\^{(o\t' ^ 0) is 
analytic in the Ihp of (O, T.^'''{(o) = E'l^(f' = 0\(o) is analytic in the uhp of (O, E^^((»)') = 
Z^l^ = 0) is analytic in the Ihp of co'. 

5. Equilibrium properties of the Green and vertex functions 

We now focus on the case that the temperatures T^. = 1/Pk and chemical potentials /x^ which 
determine the reservoir density matrices p/es' in (|9]l are all equal, 

n^T, iik^ii, k^l,...,M- (69) 

In the limits of infinite reservoir size and infinite transient time, the system is then described 
by a global equilibrium distribution: 
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^(-i)"TrpH(f=0)rc4'i''(fi)---4(''^(f'i) for V« » -'o ^ - (70) 
with the static equilibrium density matrix 

PH{t = 0)=pH = —^ , Z = Tre-'^(«-''^). (71) 

Here H is the total Hamiltonian ^ including the quantum system, reservoirs and coupling, 
and is the total particle number, which is conserved, [H,N] — 0. A proof of ( iTOl i for 
coordinate space Green functions has been given recently in ll20l for the particular case of 
a spin impurity model. 

In the present section we imply (|70] i to hold and can therefore replace the expectation 
value (fl4l i used in ( fTOl i and dTO by 

{A{t)) ^Tr pHA{t), (72) 

where the Heisenberg picture now has the reference time f = 0. 

We just described how a thermal equilibrium situation is generated by the coupling to 
large equilibrated reservoirs. In this way, equilibrium and the stationary nonequilibrium 
introduced in section |2] are embedded in a common framework. However, the following 
discussion holds also for equilibrium interacting bulk systems. The only preconditions are 
the equations (fTTT i and (|72] i, where H is the total Hamiltonian. 



5.7. Kubo-Martin-Schwinger conditions 

Any operator A (f) in the Heisenberg picture fulfils 

A(f - ij3) = el^"Ait)e-P" = p^j^t^^^ A{t)t-P>'^ pn- (73) 

Using the cyclic property of the trace, we can reshuffle two operators at times f^, tg in a 
correlation function, 

TxpHA{tA - i^)B{tB) = TrpHB(fB)e^''^A(f^)e-^''^. (74) 

We are going to apply this relation, also known as Kubo-Martin-Schwinger (KMS) 
condition lfT4l [TSll , to the n-particle Green function. Given a set of contour indices j = 
{ji,...,j„) and times t = (fi,...,f„) define f — ij3+(y') as the set of times obtained from f 
by adding (— ij3) to every fj. with j/^ ~ +. For example, given n —4 and j — (—,-!-,+,—), 
we have t — ij3+ (j) = (fi , f2 — ij3 , f3 — ij3 , f4). Let us prescribe that Tc orders operators with 
complex time arguments according to the real parts of their time arguments. Identifying A 
in (|74| | with the block of anti-time ordered operators from the + -branch and B with the block 
of time ordered operators from the — branch, we find 

G^l^(f-i^+(;)|f'-ii3+(/)) = i;'''''eP^'^''&jl,(t\t') ^ r'''^^''""'&j^,m (75) 

with G being defined in (l4TT l and m^l-'' denoting the excess of creation over annihilation 
operators on the +-branch, 

^ 1- 1- (76) 

k=l n: k=l n: 

j'k=+ 



fl(^'' ) (f 1 ) . . . fl(^") (f„ )fl(^"' {t'„)... fl(^"' * (f ; ) (77) 



The factor 11,^' ' in (iTST i is due to the different parity of the permutations which are required 
to sort 
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according to and to sort 



(78) 



according to 7^. This difference can be computed from the number M^l-' of pairwise swaps 
of adjacent operators needed to exchange the ordered block of operators from the +-branch 
with the ordered block of operators from the — branch. 



In ( fTsT l we have used that 



L 1+ L 1 



(79) 



(80) 



which is a consequence of having an even number 2n of indices. The rewriting in terms of 
m^l-' turns out to be more convenient in the following computation. 

When transforming ( fTSl l to Fourier space, we substitute on the left hand side — ij3 — > t^, 
— ij3 — > f{ for all time arguments on the +-branch and find 



,/3AJlJ (ffljci)') 



o-i/3 



n 


n d4 


n 

.k-jk=- 


n 

k:j[ = - 



q\q'\ I > 



G'\\{co\co') 



where /^^^'' {o}\o}') is defined as 
{(o\(o') = 



k-jk=- 



(81) 



(82) 



k-4=+ 

As next step we shift the path of integration for all time arguments on the +-branch by (+ij3 ) 
on the left hand side of dSTl i in order to obtain a standard Fourier integral along the real axis. 

This simultaneous shift of integration paths actually does not change the value of the 
integral, as we explain now. In principle the shift of a single time argument will change the 
value of the integral because G^'^ (f |f') is in general not analytic in any single fj. or f[; branch- 
cuts occur due to changed time ordering whenever a creation and annihilation operator of 
the same state and on the same branch of the contour have coinciding real parts of their time 
arguments. Suppose e.g. there exist ki,k2 G {I, . . . ,n} with ji^ — ff,^ — + and qf;^ =q'i^^ =: ^o^ 
that is an annihilation and a creation operator of the same state qo are both situated on the +- 
branch. For real times f/ii jf^^ '^i^^ hi approaching from above, tj^^ = the contour 

ordered operator product appearing in ( fTOl i will have the structure 



(4Ko(4)---, 



where the dots represent the other operators and a prefactor ^ if necessary, 
approaching from below, t/^^ ~ t[^ ~ 0+, the structure is 

c . . . fl^o (i (4 ) . . . = c ■ • • [1 + c«Jo K )«'?o (4 )] • • • • 

Here we use the commutation relation ( fT2] i. The results ( |83] ) and ( |84] | differ by C, so that the 
Green function is discontinuous in . When complex time arguments are allowed for, the 
discontinuity is stretched to a branch-cut in the complex f^-j -plane as sketched in figure |3] 



(83) 
But for f^-i 

(84) 
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Figure 3. If an annihilation operator at time tt^ and a creation operator at time t'j^^ are of the 
same state and are situated on the same branch of the contour, then the Green function G(t\t') 
has a branch-cut along Re(?t| ) = Re(?j.^). In this case the Fourier integrals of G(t\t') over t^^ 
along the paths Vi and are not equal. 



However, we do not only shift the integration path of a single time argument but those 
of all time arguments on the +-branch simultaneously. If the dependence on the n+ time 
arguments on the +-branch is expressed through a centre time 



1 

tc = — 

n+ 



(85) 



and (n+ — 1) relative times, then G is analytic in fc^ a simultaneous change of all times 
on the +-branch by the same value does not affect time ordering. Hence by shifting the 
integration path of tc while keeping all relative times unchanged we do not alter the value 
of the integral. (Note that the dependence on the relative time arguments might remain 
nonanalytic.) Therefore equation dSTT l is equivalent to 

g/3A^I/(«|«')G^I/(e,|o)') = («!«')• (86) 

5.2. Fluctuation dissipation theorem (FDT) 



Equation ( 186b formulates a relation between G and G. Since G is a rather artificial object of 
no direct practical use, it is desirable to eliminate it from (l86b in favour of G. 
In the one-particle case, we can make use of (l43b and obtain 

g/3A^1/(«,|«,')G^I|^,(a)|fi)') = r''' G^J^,{(o\co'), n = l. (87) 

Plugging in j,/ = =F yields four equations, two of them being tautologies and the other two 
being equivalent to 

G<{(0)^i;t-f^^"'->'^G>{(0). (88) 

The Green functions depend only on one frequency argument due to time translational 
invariance. A transformation to the Keldysh basis yields the well-known FDT, which 
formulates a relation between the correlation function G^ and the response functions G^^^ 
and G^\ 

G'^(a)) - [1 + C2«f (« - PL)] [G^'\(o) - G^'io))] . (89) 

Here 

«?(«)-^^ (90) 
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is the Bose (i^ — +1) or Fermi = — 1) function, and 

fcoth^, ifC = +l, 
1 + C2«^(«) = <^ I (91) 

[tanh^, .fC^-1. 
The FDT for the setf-energy can be obtained by combining ( |89] | with dJTl i - ( |39] ). 

E'^(«)= [1 + C2«^(a)-M)] [e'^^'M-E^^M] . (92) 
Transforming this resuh to the contour basis leads to an analogue of dSSl l. 

E<(a)) = Ce"'^('"-'''E>(«). (93) 

5.3. Time reversal 

A direct generalization of (l43T l to the multi-particle case does not exist. In order to eliminate 
G from ( |86] l. we should resort to time reversal. 

The time reversal operator is antiunitary, that is (\//|©0) = (0|0^y/) and 00^ = 
0+0 = 1; for details see e.g. |,21|. Its action either on a single-particle state \q) ~ \x,m) 
in the eigenbasis of position and s^ component of spin or on a state |p, m) in the eigenbasis of 
momentum and s- is given by 

0|a;,m) =e'''">,-m), (94) 
@\p,m) ^t'''"'\-p,-m). (95) 

Since the vacuum state |vac) is time reversal invariant, 0|vac) — |vac), the transformation of 
creation and annihilation operators turns out to be 

0flj0+ = fl^^, 0a^0+ = fl0^. (96) 
Explicitly, in the position or momentum basis it reads 

0a+ ,^0+ = &''""ai_„^, 0flx.m0^ = e-'™fl.,-,„, (97) 

0fl;„,0+ = e''^'"fl+_p @ap,,n@' = e-"^"'fl_p,_,„. (98) 

We will designate a time reversed state q or operator A by a tilde, 

q = @q, A = 0A0^ (99) 

If {s} is a complete set of orthonormal many-body states, so is {0+5}. (Completeness of 
{0+s} follows from 

(V/| [^|0M(0^-^ll m =Y.{s\©w){®<^\s) = (0</)|0r> = (rl</>>, (100) 
which holds for any i^, 0.) Hence, for any operator A 

TrA = ^(0+s|A|0+s) =^(.9|0A0+|.9)* = (TrA^*. (101) 

The time reversed density matrix is given by 

p^ = 0le-/^(«-M^)0+ = ie-/^(^-''^) (102) 

with 

Z = Tre-'5(«-^'v) = (Tr0e-^(^-''^)0+)* -Tre-/^(^-''^), (103) 
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where we used N ~ N. Hence 

Applying time reversal to an operator in the Heisenberg picture yields 

@A{t)@^ ^ ©e'^'Ae-'^'ef ^ e-'"'Ae^' = A(-f)|^, 



(104) 



(105) 



where the subscript \^ indicates that time evolution is induced by H on the right hand side. 

From ( llOll l. ( 1104b and ( I105I I we conclude for the «-particle tilde Green function defined 
in that 



]...ay (fj) 



L ''1 



= [-ir\Trpf^TjJl\-t,)\f^ 



H'i)l5 



(106) 



Here q = {qi, . . . ,q„), and in the last two lines the subscript \^ indicates that both the time 
evolution and the equilibrium density matrix are determined by H. For the very last step, we 
use the rule for complex conjugation formulated in (|47] |. In Fourier space, equation (IIO6I 1 
reads 



(107) 



5.4. Combining KMS conditions and time reversal 

UnUke the property (|43] l of G which is restricted to the one-particle case, equation jlQlj holds 
for any n and provides a general possibility to eliminate G from (|86] |. Using it we find 

^/3A.I/(H«>')GA/(e,|co') = C"'' G^^\l{co'\co)\f^, (108) 

which is a generalization of (|87] | to arbitrary particle number n. We recover formula dSTb 
from dTOSl l by means of 

Gjj{,(a)|a)') = G^l|l(«'|a))|^, « = 1, (109) 

which follows from combining the Fourier transform of (l43T l with (I107l i. 

If one transforms equation (1108b to the Keldysh basis, the result can be considered 
as a multi-particle FDT since the components of the multi-particle Green functions in the 
Keldysh basis have a distinct interpretation in terms of response and correlation functions fT9l. 
We perform this transformation in section \5M restricting ourselves to a class of physically 
interesting systems which allow for certain simplifications of the multi-particle FDT. 

Before deriving an n-particle identity for the vertex function analogous to dlOSb . we note 
that the bare two-particle interaction vertex given by (fTSl l satisfies 

^/3a/U(«,»,/I;(^,|^) = r'''v§{co\co%. (110) 

In the case ^ {( | )5 (+ + I + +)} this equation is automatically fulfilled because 

both quantities equal zero. For {f\j) G {( | )5 (+ + I + +)} equation (II 101 ) states 

V^^)^^{C0'\C0)=vl^l,{(0\(0% if /i=;^ = yi=;2 = T, (111) 
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which follows from 

{q'W2H\<im) + l^\q2qi)) = ((^W^l + C(^k'i I)0^0v0t0|^i^2) 

= {qiqi\v{m2) + Q^2^i)). (112) 

where v = 0v0^ 

With help of jl lOl i we can now prove the identity 



(113) 



by identifying one by one the diagrams contributing to each side of it. The values of the 

individual diagrams are given by The prefactors e^^' and C"'' ' appearing in dl 131) 

can be split and distributed among the noninteracting single-particle Green functions and bare 
two-particle vertices in ( l20l i by making use of the properties 

A>'l^'(a)'|a)) = A^''l'(a)'|v)+A'l^'(v|a)), (114) 

m-''!' + m'l-'', (115) 

which can be easily generalized to mixed particle numbers and chains of arbitrary length, like 
e.g. 

„/l72l)i;2 ^ ^.7'l>2l'l'2 +m'2l'2 +m'l'2l-''l-'2. (116) 

According to ( 1108b and ( II lOl l the constituents of both sides of (II 131 ) can then be readily 
identified. 



5.5. Systems with special behaviour under time reversal 

For a further evaluation of ( 11081 ) and (1113b . we should establish how the system behaves under 
time reversal. We restrict the following discussion to the class of systems which satisfy 

for all indices and time arguments. Note that equation (II 17b is not only a condition on the 
Hamiltonian but also on the basis of single-particle states. 

We mention some examples which demonstrate the extent of this class of systems. A 
trivial example for the validity of (11 17b is given when the Hamiltonian and the single-particle 
states are time reversal invariant. 

As a less trivial example we consider the single impurity Anderson model (SIAM) lfT6l 
in a magnetic field B. This model consists of an electronic impurity with on-site interaction 
which is coupled to a conductor. The corresponding Hamiltonian reads 

HB^Y^{£Q + <yB)a'„ac + Y^ j d^ p Epc'^p^Cp^a 
a a •' 

+ i d^p [Vpalcp^a + VL-c] +f/a|flt«|fl|. (118) 
a •' 

Here CT = ±5 =tii denote the eigenvalues of the single-particle spin i.., and c\c are creators 
and annihilators of states in the conductor. 

As a single-particle basis to be used in (II 17b we choose the common eigenbasis of 
position and s~, \q) = |a;, o), where the values of x are considered to label the states of both 
the conductor and the impurity. 
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In order to verify (11 17| i for this Hamiltonian and this choice of basis, we make use of 
the unitary transformation D. — e""^'^' rotating spin space by n around the x-axis, being the 
;ic-component of the total spin, 

^\x,a) = -i\x,-a). (119) 
Using ( |98] l and the fact that Sp and Vp do not depend on the sign of p in (11181 ). we find 

Hb = H^B = ^HbQ.^ . (120) 
Since a combination of i94[ and ( II 191 ) yields additionally 

=0|a;,(T) =e"'(^+2)i2|a;,(7) =e'''('^+i)i2|^) (121) 

we conclude 

with £lq — i^qi,- ■ ■ ,^qn)- In (1122b we exploit that the Green function is nonvanishing only 
for 

N N 

L - L (123) 

k=l k=l 

due to spin conservation. It is straightforward to derive from ( fTOl i. ( fTTT i and ( |72] | that 

Gng\nq'\nHat ^^iW- (124) 
[Note that ( 1124b holds in general for any unitary operator £2.] Combining (1122b and ( 11241) 
we conclude that the condition ( II 17b is fulfilled for the impurity model dl 18b and the single- 
particle basis 1^) = \x,a). 

On the full analogy it can be shown that dl 17b is also valid for this impurity model and 
the single-particle basis \q) — \p, a). In ( 1120b - ( 1124b . £2 simply has to be replaced by £211, 
where 11 is the parity transformation, 

n|p,a) = |-p,a). (125) 

As an example of a system which does not fulfil (II 17b we consider the Hamiltonian 

HA = jA\al^ (126) 

for spinless particles with charge e moving in the field of a vector potential A{x). Let us 
choose the position eigenstates as the single-particle basis. Since 

Ha=H_a, (127) 
we conclude that for a generic vector potential 



^ = G 
Ha 



7^ G I ' / 



(128) 

Ha 



Obviously, equation (II 17b does not hold for this choice of Hamiltonian and single-particle 
basis. 
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5.6. Generalized fluctuation dissipation theorem for systems satisfying equation jllTIl 
When relation jl 17b holds, equation ( 11081 ) takes the form 

g/3A^I/(«,|«,')G^1/(e,|co') = C"""Gf l|;(a)'|a)). (129) 

Additionally, equation (11 17t can be expanded in orders of the two-particle interaction. From 
the first-order term of this expansion of the two-particle Green function, we infer that the bare 
two-particle interaction vertex complies with 

v%ico'\co) = d!^{co'\co)\^. (130) 
When combined with ( II 10b it brings about the identity 

g/3A/U(«,»-/U(^/|^) - C'"''''v^'lJ(«l«')- (131) 
Similarly to the derivation of dl 13b we deduce from (1 129b and ( 1131b that 

g/3A/U(«,»y;j;(^,|^) = C"'''''^K(«I«')- (132) 

v I w '■frf 

It is convenient to combine ( |129b and ( 1132b with ( l48T l and (|50] |, respectively, in order to obtain 
the equations 

g/3A^I/(«|<«')G^I/(e,|e,') = (-i)«C'"''''GjJ(to|to')*, (133) 
^/3a/u««»^;|.(^,|^) = (134) 

where state indices and frequency arguments are now identical on both sides. 

We are going to transform ( 1133b and (1134b to the Keldysh basis in order to achieve a 
multi-particle FDT. In what follows the state indices are omitted for brevity. According to 

G^l^" = ^ z)^1«G"l"'(D-i)"'l^" (135) 
a.a' 

every component of G^'-' is expressed as a certain linear combination of components of G"'"'. 
Let us define 

Gi^j' ^ £ I !)«'!/ (136) 

a,a' 

for e = ±1 (cf. 112113). It is obvious that 

Gj\j' ^G^y +G^1^' . (137) 
Using the relations 

dJ\" = (_i)i:*«AZ)/l«^ (138) 

^ (_l)E*«A-(£)-l)«l/ (139) 

which follow from ( |29] l - (l3Tb . we find additionally that 

GJ\j' ^ G^y - G^l^' . (140) 

Plugging ( 1137b and ( 1140b into ( 1133b and splitting the latter into the real and the imaginary 
parts we obtain 

sinh ^^^'^'j'"!'"'' ReG^'l-'',,(a)|a)') - -cosh ^^^''V"'' ReG^''^" .,,(«|«'), (141) 
cosh 2^^^^ ,.,{C0\(0') = -sinh^^^^^ ImG^I^' ,,(«|«'), (142) 
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where 



e,;'^' = . (143) 

Utilizing ( |9l1 l we finally achieve the representation 

RtG'^^l,{(o\(o') = -[l + C2nf (A^'l^"((o|t(j'))j ReG^''-'"^.|.,(a)|«'), (144) 

lmG^''^"^.|^, {(o\(o') = - [l + C2n^(A^''^''(co|a)'))j ImG'''^^, (to|to'). (145) 

These two equations express relations between certain linear combinations of the components 
of G"l" . Since the components of G"\" can be understood in terms of higher order response 
and correlation functions |fT9l . equations (1144b and ( |145t constitute a multi-particle FDT. 
Both equations can be evaluated for 2^" different realizations of the indices {j\f)- However, 
the underlying equation ( 1133b is equivalent for G^'^ and for G-''-'' as can be inferred from 
fj^Jlf = and Aj\j' {co\co') = -A^^j' {co\(X)'). Therefore only half of the 2^" equations ( 11441 ) 

are independent. The same applies to (1145b . 

For bosons the function nf=+ (A-'l-''(a)|a)')) diverges in the special case 
(which occurs for all contour indices being — or all being +), and one then should rather 
use (fTTTT l and (fT42l l with cosh[j3A-'l-''(a)|«')/2] = 1 and sinh[j3A^I-''(«|a)')/2] = instead 
of dm and (fT45] l. 

The transformation of (11341 ) to the Keldysh basis is completely analogous to that of (1133b . 
Defining for e = ± 1 



(146) 



a ,a 



we end up with 



^1 

Im/l^.,^/«»: 



l + C2«^(A^"l-''(a)'|a)))jRe/l^.,|/a)», (147) 

^1 

1 + C2n^{Aj'\j{co'\co))] lm/j(.i(o'\(o), (148) 



where 



l+mJ y 



(149) 



[cf. ( Il43b 1. Like in the case of Green functions, each of equations ( |147b and ( 11481 ) contains 
only 2^"^' independent relations. 

Let us illustrate equations ( 1144b . ( 1145b and ( |147b . ( 1148b by some examples. In the one- 
particle case, n = 1, we evaluate (1144b and (1145b for {j\f) — {—\ — ) and {j\f) = {—\+)- Since 



i(Gi|2- 





f G2|2) and G-I+ = 




2 2 


g: 


= i(Gi|2 + G2|') = 


g;I+ 


2 2 


g:I+ 


= i(Gi|2_G2|') = 



i(G^|2- 



■^211 



G I ), we have 



(150) 
(151) 
(152) 
(153) 
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Furthermore, we use 

^-\-{(o\co')^0, er'" = -C, (154) 

A-l+(a)|«') = «'-Ai, er'^ = -l, (155) 

and the representation of the one-particle Green function G{(i)\(i)') = 2n5{(0~ (o')G{(o). 
Then equations (11441 ) and (I145I I evaluated for {j\f) = (-|-), (-|+) read 

ReG^{a))=0, (156) 

Re [G'^^(to)-G'^"(a))] = - [l + C2«^(« - ju)] ReG'^(a)) = 0, (157) 

lm[G^''{(o) + G^^'ico)]=Q, (158) 

ImG^ico) = -[! + C2ni;i(0 - n)] Im [G^'{co) ~ G'^'^co)] . (159) 

These formulae obviously reproduce ( l54l ) and ( |89] ) under the additional constraint G^|^/ (o) = 
Gq'\q{co) which follows from the special one-particle property ( 11091 ) and the assumption ( II 17l l. 

Evaluating ( 1147b and ( 1148b for the one-particle vertex (the self-energy), leads to a result 
similar to ( 11561 ) - ( |159b : one only needs to replace G therein by the self-energy E. 

For n = 2 the complete set of relations for the Green functions is obtained by 
evaluating (1144b and ( 1145b for eight independent choices for (./i j2|./i 72) from 16 existing 
possibilities. Let us for example discuss the three particular combinations {jij2\j'ij'2) = 
+ l --),(--!++)■ Since 



A""l""(a)i,a)2|a)J,a)2) =0, 


£2 =c, 


(160) 


A"+l""(a)i,a>2|«J,to^) = -(02 + ^1, 


e2"+'"" = +l, 


(161) 


A""l++(a)i,a>2|«J,o>^) = + «2-2^, 




(162) 


equations (11441) and (|145|) read 






ReGI"' (0)1,0)2 oj, 0)2) =0, 




(163) 


ImG!|!"' ( 0)1, o>2 0)1,0)^) =0, 




(164) 


ReG;^'""(0)i,0)2|0);,0)2) = -[1 + C2nf(-0)2 + Ai)] ReGI^'"" 


"(0)i,0)2|0);,0)2), 


(165) 


ImGI^'""(0)i,0)2|0)I,0)^) -[l + C2n^(-o>2 + Al)] ImG^^'"" 


(0)1, 0)2 1 0)1,0)^), 


(166) 


ReG^"'++(o)i,0)2|o)J,0)2) = -[1 + C2nf(o)| +0)^-2Ai)]ReG; 


:J++(0)i,0)2|0);,0)^),(167) 



ImG_j'++(o)i,0)2|o)I,0)^) = -[l + C2nf(o)( +0)2-2^)] ImG^ '++(o)i,0)2|o)J,0)^). (168) 

The functions (j^^^liiii appearing in these equations can be found by explicit evaluation 
of ( fT36] l with help of (|29| - (EB- This results in 

G;-l- = i(G'i|i'+G''|22 + G'2|i2 + G'2|2i+G2i|i2 + G2i|2i+G22|ii), (169) 

^--1-^ 1^^11|12^^11|21 ^^12111 ^^21111 ^^22121 ^^22|12^^21|22^gl2|22)^ (170) 
^- + 1-^ 1/ ^11|11_^11|22^^12|12^^12|21_^21|12_^21|21^^22|11N (^^^^ 



1 



Gr'-- = -(-Gll|12-Gll|21+Gl2|ll-G21|ll+G22|21+G22|12_G21|22 + Gl2|22)^ (172) 
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^--1++^ 1 ^^11111 ^^ll|22_^12|12_^12|21_^21|12_^21|21 ^^22111)^ (173) 
^--1++^ 1(_^11|12_^11|21 ^^12111 ^^21111 _^22|21_^22|12^^21|22^^12|22)^ (174) 

6. Conservation of the properties in approximations 

In practice the Green and vertex functions can often be computed only in approximations. 
In this section, we address the question whether such approximations are consistent with 
the exact relations presented in this review. Since the answer depends on the details of 
the corresponding method, we have to restrict ourselves to specific approximation schemes. 
First we discuss the large class of diagrammatic approximations, and then we analyse 
approximations within the framework of the fRG. 



6.1. Diagrammatic approximations 

Many of the approximations used to compute the Green or vertex functions are given by taking 
into account only a subset of all diagrams contributing to that function. This may be either a 
finite set like in perturbation theory or an infinite one like in the random phase approximation. 
We call this approach a diagrammatic approximation and assume that apart from choosing 
only a subset of diagrams no further approximation is done. Since all discussed relations 
are linear in the vertex functions it follows that if two distinct diagrammatic approximations 
satisfy a given relation each, so does their sum. 

Consider a diagrammatic approximation to a multi-particle Green or vertex function. 
The condition for the approximation to fulfil the relation ( l45l l or (l46l l under permutations of 
particles is that the set of diagrams taken into account is invariant under those permutations. 
For instance neither of the diagrams in figure |4] alone fulfils the relation (|46] |. The 
approximation given by the sum of diagram (a) and (b) in figure |4] satisfies a part of ( |46] l, 
namely 

7;;;|^'(i^a)>)=r^:j;(«'i«). 

The same is true for the sum of (c) and (d). The sum of (a) and (c) satisfies 

^:|;^(«'|Pa))^r^i(«'|a)), (176) 

so does the sum of (b) and (d). The sum of (a) and (d) satisfies 

Yll\li^{P(o'\P(o) = Y!,^{(o'\(o), (177) 

as does the sum of (b) and (c). Only the sum of all four diagrams fulfils the complete 
relation (|46] |. In practice one will exploit the symmetry imposed by ( |46] | to reduce the number 
of components which have to be evaluated directly. For instance, it would be sufficient to 
calculate one of the four diagrams in figure H] and use the result to derive the other three 
diagrams via ( |46] |. 

Since complex conjugation of a Green function interchanges creation and annihilation 
operators it reverses the direction of lines in diagrams. This results in the following criterion 
for a diagrammatic approximation to comply with equation ( l48T l or ( fSOb [or equivalently ( l52b 
or (l53Tl1: the set of diagrams taken into account has to be invariant under reversion of the 
directed lines combined with the simultaneous exchange of the external indices j with /, q 
with q' , and o) with (£)' . For example, reversing the direction of all lines in figure 2] (a) and 
interchanging 1 with 1' and 2 with 2' reproduces the same diagram. Hence, an approximation 
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(a) (b) (c) (d) 




Figure 4. Examples of diagrams contributing to the two-particle vertex function 7i'2'|i2- The 
external indices represent contour index, state and frequency, e.g. 1 = (j\,q\,0}\). 

to the two-particle vertex function consisting only of this single diagram complies with ( fSOl l. 
The same holds for figure|4](d). The diagram in figure|4](b) is mapped onto that of figure|4](c) 
by complex conjugation and vice versa. The approximation given by their sum fulfils ( fSOl l. 
Similarly to the rules for permutation of particles, one will use (|48] | or ( ISOl l in practice in order 
to reduce the number of components one has to evaluate directly. 

The theorem of causality for the vertex function has been proven in section |43] for each 
single diagram; the theorem of analyticity is a direct consequence of it. Therefore, both 
theorems hold in any diagrammatic approximation to the vertex functions. Given the causal 
features of the noninteracting single-particle Green function, the theorem of causality for the 
interacting multi-particle Green function can also be proven to hold for any single diagram. 
The proof is analogous to the one for the diagrams of vertex functions given in section 14.31 
That's why the theorems of causality and analyticity for the Green functions hold in any 
diagrammatic approximation as well. 

The proof of the KMS condition ( II 13l l for the vertex functions has been performed 
diagram per diagram. The same applies to the special form ( 1132b which brings about (I147l i 
and (I148I I. Hence ( 1113b and (1147b . ( 1148b hold in any diagrammatic approximation to the 
vertex functions. Again, an analogous proof for the diagrams of the interacting Green function 
can be formulated, showing that equations (1108b . (1144b . and ( 1145b hold in any diagrammatic 
approximation to the Green functions. 



6.2. Functional renormalization group 

While diagrammatic approximations automatically preserve causality and the KMS 
conditions, these features can be much more difficult to maintain in other approximation 
schemes. We briefly sketch complications for the fRG approach which is based on introducing 
a flow parameter A into the free single-particle Green functions; typically this is done in such 
a way that degrees of freedom at energies below A are suppressed in order to regularize low 
energy divergencies which may arise in perturbation theory. The vertex functions, expressed 
as functionals of the free single-particle Green functions, depend on A and can be computed 
as solutions of a coupled set of flow equations Il22l |23l 1241 . 

Let us consider the flow of the one-particle irreducible vertex functions 7,, a governed by 
the hierarchy of the fRG equations 

d v'lv v'^'IV'^r" /i-70\ 

^YlX =C72;A 5a;A|A', (178) 

d v;v^|VlV2 _ ^ VJV^A'|VIV2A , v;v;|Ai;i2 ^^l'^|VlV2 

^72;A -hYi-A ■^A;;i|A' + 72;A WiA, |A,"^A;A2|A^ 72;A 

+ \^A-li\AJ,^A;l2W l^yi-A 72;A + 72;A 72;A J + (,Oa O ^Aj | , (1 

(180) 
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where 71 a = Ea, and the muhi-indices v and A comprise frequency (time) argument as well as 
the state and the Keldysh indices; a summation convention applies to multi-indices appearing 

twice in a product. The dots in (II8OI 1 represent the flow equations for 73;a,74;a, The 

notation S'a stands for the so-called single scale propagator which is defined by 



The equations (I178I I - (II8OI 1 can be systematically obtained from the vertex generating 
functional explicitly depending of the flow parameter A [^T TOl . An equivalent derivation 
is based on a purely diagrammatic approach 1 1 1 1 in which A is explicitly introduced into the 
free propagators G^"^. 

Since the flow of 7„;a is determined by 7i;a, ■ ■ ■ , 7n;A and 7„+i;a, the exact flow 
equations ( I178l4180"l l form an infinite coupled hierarchy. The final solution of the full hierarchy 
is exact and therefore satisfies all relations discussed in the previous sections. However, 
in nontrivial practical problems one cannot solve the full hierarchy but has to resort to 
approximations. Typically one truncates the hierarchy of flow equations which means setting 
d7«;A/dA = for some n, such that the flow equations for 7i;a, ■ ■ ■ , 7n-i;A form a finite closed 
set. Then it is no longer guaranteed that the final solution of the truncated flow equations 
satisfies the exact relations in question. It turns out that the choice of the flow parameter has a 
decisive influence on their conservation. While arbitrary choices of the flow parameter do not 
affect the properties under permutation of particles and complex conjugation, special care is 
required in order to maintain the theorem of causality and the KMS conditions. 

Assuming that the only approximation done to the flow equation is a truncation of the 
hierarchy we find the following sufficient conditions for the construction of conserving flow 
parameters. (1) If the free propagator Ga ' furnished with the flow parameter A satisfies the 
theorem of causality, then the vertex functions resulting from a truncated set of flow equations 
fulfil the theorem of causahty throughout the flow and in the final solution. (2) If Ga ^ satisfies 
the fluctuation dissipation theorem ( |89] l (with the factor [l + ^2n^{co — /i)] not depending on 
A) then the vertex functions resulting from a truncated set of flow equations fulfil the KMS 
conditions throughout the flow and in the final solution. 

For the proof of these rules we recall that the theorem of causality and the KMS 
conditions for the vertex functions have been derived in sections |431 and 15341 for each single 
contributing diagram on the basis of the corresponding property of the free propagator G'*'' . 
In reference 1 1 1 1 the flow equations for the vertex functions have been derived from the 
individual diagrams contributing to d7„;A/ dA, which contain lines Gj^' and one line dGj^^ / dA. 

If Ga^ and dGA satisfy the theorem of causality, then the same arguments as applied in 
section l43] can be used to prove the theorem of causality for each single diagram contributing 
to d7„;A/dA. Truncating flow equations means taking into account only a subset of these 
diagrams and it thus has no influence on the conservation of causality during the flow. Note 
that the theorem of causahty for dGA follows from that of Ga' since the derivative is a 
linear operation. Concerning the KMS conditions we can argue analogously. 

As examples for these considerations we discuss different flow parameters used in the 
literature and propose also a new one. In reference ifTOl the flow parameter is introduced into 
the free Green function by setting 



Since the step function &{z) has a branch-cut along the imaginary axis of z, the analytic 
properties of G^^-{o}) are in contradiction to the theorem of causality; G^^\(o) being 




Gf («) = 0(|«|-A)gW(«). 



(182) 
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nonanalytic in the upper half plane implies that G^^^^^{t\t') — G^'''^''(f |f') does not vanish 
for t' > t. As a consequence the vertex functions constructed from diagrams with lines given 
by of^ instead of G^''^ do not exhibit the correct causal features. When the exact RG flow 

stops at A = 0, causality is restored due to cf^^^ — G^^\ However, for nontrivial models 
the flow equations for the vertex functions can be solved only approximately. Restoring 
causality for approximate solutions of the flow equations was encountered as a major difficulty 
in reference [10|. 

This problem can be solved if one chooses a flow parameter such that Gj^'(a)) has 
the same analytic features as G^^\(o). Then the causal features of the vertex functions 
are preserved during all the RG flow, even for truncated flow equations. An example of 
such a flow parameter is the imaginary frequency cut-off proposed in reference lillj which 
manipulates the particle distribution functions of the reservoirs by 

where a)„ = 2nnT or O),, ~ {2n + l)nT denote the Matsubara frequencies for bosons or 
fermions, respectively. This cut-off affects only the Keldysh component of the free single- 
particle Green function, but not the retarded or advanced ones which contain only dynamical 
but no statistical information, Gf-^^^^^''{co) = G''^^^^^^ {co) . Hence, the causal properties of 
free single-particle Green function are unchanged and causal features of the vertex functions 
are conserved. However, since the A-dependent distribution function in ( |183t is not the 
thermal one, the single-particle Green function G^^^ does not fulfil the fluctuation dissipation 
theorem for A 7^ 0. As a consequence the KMS conditions for the vertex functions are 
not satisfied until the exact flow terminates at A = 0. Solvable approximations to the flow 
equations will in general not restore the KMS -conditions even at A = 0. (The simple static 
approximation used in reference fTT) yet satisfies the KMS -conditions in a trivial way.) 

Examples of flow parameters which violate neither causality nor the KMS conditions are 
given by the momentum cut-off and by the hybridization of the system and reservoirs. 

The momentum cut-off is introduced for extended systems by defining 

^(WAv/K^^^ ^ 0(|e,| - A)Gf ^^'/^^/'^(o)), (184) 

where the wave number k labels a single-particle eigenstate with eigenenergy Obviously, 
this cut-off does not affect the frequency dependence of G'"'; hence, it preserves causality and 
the KMS conditions. 

We propose to use hybridization as a flow parameter which is useful for local systems 
coupled to external reservoirs. After integrating out reservoir states, the Green function of 
a local system acquires a reservoir broadening F. For example, the reservoir-dressed Green 
function of the SIAM defined in (III8I1 reads 



G^^^V) ^ (185) 

where F ~ 2nY,p8{^ — ep)\Vp'^. The idea of the hybridization flow parameter consists in 
introducing an additional - fictitious - reservoir, which is assumed to be in thermodynamical 
equilibrium with the physical one, i.e. is characterized by the same temperature T and 
chemical potential fX . The flow parameter A is then associated with an additional broadening 
due to a coupling to the fictitious reservoir, so that 

G^^l'\(0) = ^ ■ , (186) 

"•^ ' «-(eo + aB) + ^(F + A)' 
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'a- A 

-,(0)K 



(CO) 



-- [l + C2n^(a)- 



■,(0)Av 



(0)) 



(187) 
(188) 



In the beginning of the flow at A = oo, the free propagator is suppressed and the vertex 
functions are identical to their bare values (except for the self-energy which acquires a final 



contribution from the Hartree-Fock diagram, E 



Ret 



;A= 



o). In the end of the flow 

(0) 



at A = the fictitious reservoir is decoupled, and we restore the original system. Since G^ in 
jl 86111 88] l satisfies the theorem of causality and the fluctuation dissipation theorem, causality 
and the KMS conditions of the vertex functions are preserved during the flow, even when the 
flow equations are truncated. A detailed study of the SIAM on basis of the hybridization flow 
is presented in [25|, where it is also generalized to nonequilibrium situations. 



7. Conclusion 



Currently the interest in nonequilibrium dynamics and in the crossover between equilibrium 
and nonequilibrium properties is enormously increasing in the field of condensed matter 
physics. Keldysh formalism provides the possibility to treat equilibrium and nonequilibrium 
quantum systems on equal footing. Numerous equilibrium techniques based on Green 
or vertex functions benefit a lot from adaptation to Keldysh formalism. For its correct 
implementation the knowledge of real time (frequency) features of the multi-component 
Green and vertex functions appears to be mandatory. The present review lists general 
properties of these functions in equilibrium and nonequilibrium stationary state. 

The relations for the Green functions have been derived directly from their definition as 
expectation value of contour ordered operator product. The corresponding equations for the 
vertex functions have been obtained from the rules for the evaluation of irreducible diagrams. 

We have surveyed the relations for permutations of particles (145 1 146 1 and complex 
conjugation (|48ll50t and cast them into the form convenient in practical use. In section |43] we 
have formulated the causality theorem and identified the components of the Green and vertex 
functions which are analytic in a certain half plane of the corresponding frequency arguments. 
These analytic features should be maintained throughout any approximation to the Green and 
vertex functions, otherwise the fundamental property of causality is violated. 

A major part of the review has been devoted to the study of the KMS conditions which 
hold in equilibrium. We have worked out the constraints which equilibrium imposes on the 
multi-particle functions. They are given by the relationships (I108II1 13l l to the multi-particle 
functions of the time reversed system. The connection to the FDT for single-particle functions 
has been extensively discussed. In particular we have determined the reason for its special 
form ( [89] l for n — I: it follows from the peculiar single-particle property d43T l. Furthermore, 
we identified an important class of systems with special behaviour jllH under time reversal 
which allow for a simpUfied formulation ( ll44lll45l[T47lll48l l of the multi-particle FDT. 

Finally we discussed how the properties in question relate to diagrammatic 
approximations and to approximations within the fRG. For diagrammatic approximations we 
found that the set of diagrams taken into account should be invariant under permutations 
within the incoming and within the outgoing external indices and invariant under reversal of 
the directed lines. Then the relations for permutation of particles and complex conjugation 
hold and reduce considerably the number of independent components which simplifies a 
treatment of a problem. Causal features and KMS conditions are preserved automatically 
in any diagrammatic approximation. However, approximations to the flow equations of the 
fRG often destroy causality and the KMS conditions. These properties can be preserved when 
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the flow parameter is appropriately chosen in such a way that it respects the corresponding 
properties of the noninteracting single-particle Green function. 
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